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Abstract 

The total energy in the two-phase chiral bag model is studied, 
including the contribution due to the bag (Casimir energy plus energy 
of the valence quarks), as well as the one coming from the Skyrmion 
in the external sector. 

A consistent determination of the parameters of the model and the 
renormalization constants in the energy is performed. 

The total energy shows an approximate independence with the bag 
radius (separation limit between the phases), in agreement with the 
Cheshire Cat Principle. 
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The hybrid chiral bag IB is an effective model to describe the behavior 
of strongly interacting baryons. In this model, color degrees of freedom 
are confined to a bounded region and coupled to a bosonic external field 
(skyrmion) through boundary conditions. 

These two phase models are intermediate between two successful descrip¬ 
tions of baryons: bag models [4| - with QCD degrees of freedom at short 
distances - and Skyrrne model i i- a, an effective (non renormalizable) 
nonlinear sigma model, useful when the low energy properties of baryons are 
considered. 

An interesting feature of Chiral Bag Models (CBM) is the appearance 
of the so called Cheshire Cat Principle (CCP) B 1, according to which 
fermionic degrees of freedom can be replaced by bosonic ones in certain re¬ 
gions of space, the resulting position of the limit of separation between the 
two phases having no physical consequences. 

In 1 +1-dimensions, the Cheshire Cat behavior follows from the bosoniza- 
tion of fermionic fields 0 . In the 3 + 1 case, topological quantities, such as 
the baryonic number, have a similar behavior but, for non topological 
ones, the CCP is expected to be only approximately valid. 

In what follows, we will study the energy of a four-dimensional hybrid 
model consisting of quarks and gluons confined to a spherical bag plus a 
truncated exterior Skyrrne field in a hedgehog configuration. It is our aim to 
study the dependence of the total energy on the size of the bag, thus testing 
the Cheshire Cat hypothesis. 


In bag models, quarks and gluons are confined to a bounded region that, 
in our case, will be taken as an static sphere of radius R. Adequate boundary 
conditions are imposed on the field so as to guarantee the vanishing of the 
flux of color current to the external sector. 

In the MIT bag model [3| the boundary conditions for fermions and gluons 


are 


B*\b = -(Z + i to* 

n^\ B = 0 . 


= 0 


B 


( 1 ) 

( 2 ) 


The Dirac operator for the fermionic field, together with the boundary con¬ 
dition ([!]) define an elliptic boundary value problem [1C], Moreover, i) gives 
a vanishing current at the boundary, T/mT = 0. 
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But the axial symmetry, present in QCD with massless quarks, is broken 
at the boundary when condition (|l|) is satished. To avoid such unwelcome 
behavior, CBM [JJ| have been proposed, which amount to imposing on the 
fermionic held the boundary condition 


AT L = 


B 2 


\\l + i fie~ ie ^ Hh5 ) T 


B 


= 0 , 


(3) 


where r* are Pauli matrices, and adding an external sector to the bag. 

Through this boundary condition, the fermionic held is connected with 
the external held represented by the “chiral angle” 8(R). As it will be shown 
later, the external phase can be described using the Skyrme model. 

To obtain the energy in the chiral bag model we will proceed in steps, 
using some results previously obtained in references 
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In the hrst place, we will introduce the difference between Casimir ener¬ 
gies of chiral and MIT fermionic bags. It corresponds to the zero temper¬ 
ature limit, T —> 0, of the results presented in [|TJ. As a second step, we 


will study the reference vacuum energy, i.e., the MIT Casimir energy, which 
is the T —> 0 limit of the cases analized in |12| for the fermionic and gauge 
helds. An external Skyrme held will then be introduced, so as to complete 
the two phase model (TPM). The caracteristic parameters in its Lagrangian, 
together with the renormalization constants in the bag energy, will be de¬ 
termined through physical considerations, thus obtaining the total energy of 
the TPM. 


Determinants of quotients of elliptic differential operators under differ¬ 
ent boundary conditions can be expressed as p-determinants of quotients of 
Forman’s operators |13], [L4], [15]. This leads, for a bounded euclidean time, 
to the study of differences of free energies of the physical system subject 
to two different boundary conditions. In such a way, in reference [|TT]] the 
difference between the free energies of an SU( 2) chirally symetric system of 
massless fermions, confined to a spherical region and subject to chiral and 
MIT boundary conditions respectively, was calculated. 

To construct Forman’s operator, which is totally defined by its action 
over functions in the kernel of i ft, a discrete basis of the space of solutions 
was considered. The differential operator and the boundary conditions are in¬ 


variant under the diagonal subgroup of SU(2) rota tion® SU(2) 


isospin > 


and leave 
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invariant the subspaces caracterized by the quantum numbers {k,j, l, m}, as¬ 
sociated with the eigenvalues of {K 2 , J 2 , L 2 , K z }, where K = L + S + I. In 
particular, k takes non-negative integer values. 

The T —» 0 limit of the above mentioned diference of free energies leads 
to 

Ae c (6) = R{E C (R,6) - E cMIT ] 

— \AtiKq sin 2 6 + 0.463 sin 4 6 + 0.023 sin 6 9 
47T L w j 


= 3 


1 oo /»oo 

— y dx 2v log (l + Ck{x) sin 2 9 + D k (x) sin 4 9 — A(k, x)) 
2n k=1 J o L V j 


1 r°° 

/ dx log 
2n J o 


J_ f e 2 0 < 9 < tt/2 1 

47T | (vr — 6) 2 7t/2 < 9 < tt j 

(x) -4+[2(l-a 2 (q ; )) 2 + ^1cos 2 e 


1 + 


x 4a 2 (a:)+(l— a?(x)) 2 cos 2 8 


^ _ ol{x) _ 1 


x l-\-a 2, {x ) 


(4) 


where the Casimir energies have been adimensionalized. The first and second 
terms in the r.h.s. of (|4|) correspond to the subspaces k > 1. The third and 
fourth ones, to the k = 0 case. We have used the definitions 


C k (x) = 


-2 


■x 


4 x 2 dl(x) + (p 2 - d? k (x)) 

(/ P 2 - d 2 k (x)) 2 4 x 2 d 2 k (x) + (p 2 - dl(xj) 
T ( P 2 - z' 2 ) 4 x 2 n d 2 k (x) - (p 2 - d 2 k (x)) 


(5) 


Dfz{x) = t 


(p 2 — d 2 (x)) 2 — (p 2 — v 2 ) 


4 x 2 d 2 k {x) + (p 2 -d 2 k (x)Y 


( 6 ) 


v = k + 1/2 p = Vx 2 + v 2 
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d k (x) = x— In I v (x) 
ax 


a(x 


= coth x 


a(x) = 2 a(x) 


x 


where I v (x) is the modified Bessel function. In (||), A (k,x) represents the 
hrst few tenm£] in the asymptotic (Debye) expansion of the other terms in 
the argument of the logarithm, required to isolate the divergent pieces in the 
Casimir energy. 

In obtaining (||), an analytic regularization of non-absolutely convergent 
series has been performed through the introduction of the factor p~ s , for 
3?(s) large enough, and then taking the finite part at s = 0 0 , giving rise 
to the hrst term in the r.h.s. This procedure leaves an arbitrary finite part 
proportional to sin 1 2 (6 ) ), which requires the introduction of the undetermined 
constant Kq. This amounts to the introduction, in the Lagrangian of the 


external Skyrrne held, of the counterterm |T6 


Kn 


IQlxR Jr= 


f d 2 xtr \L a L a — (n a L a ) 2 ] = ^-sen 2 d, 

Jr=R 1 > R 


(7) 


where L a = Wd a U = e~ ie ^ T ' r ^d a e ie ^ T ' r \ 

The second and fourth terms of ([|) require numerical calculations. In the 
hrst case, the sum over the index k has been cutted when the tail of the series 
becomes negligible. Integrations in the x variable have been numerically 
solved in both cases. 


As it was said before, the reference (MIT) Casimir energy for fermions 
and that corresponding to gluons have been studied in a complementary way. 
In reference |12[, the free energy for a fermionic and an abelian gauge held 
(enough for the 1-loop description of the free energy for gluons) was studied. 
Such evaluations have been made using analytically regularized traces which 
involve the Green functions of the boundary problems considered. 

The T —> 0 limit of the results in JT^| has a simple structure. As it can be 
understood by dimensional analysis, the Casimir energy (once singularities 


1 Note that, in order to isolate divergences, it is enough to retain the first three terms 

in the asymptotic expansion 0 - However, in the present calculation, we have retained 
the first six terms in the Debye expansion for computational convenience. 
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have been removed by the renormalization of the zero-loop energy) takes the 
form 

e c,M/T = RE cM it = Emit , (8) 

where Km it is the arbitrary finite part left by the renormalization procedure. 

We thus have now the total Casimir energy of the bag, including the 
correction due to the interaction of fermions with the external Skyrme field, 
represented by the chiral angle 9, 

e c (9) = A e c (6) + e c ,MiT , (9) 

up to the knowledge of the constants Kq and Km it■ As we will show later, 
they can be determined imposing physical conditions in the framework of the 
two phase chiral bag model. 


A further contribution has to be included if the total inner bag energy 
(energy of the defect, in the remaining of the paper) is studied. 

In reference ||, it has been proved that the valence quark contribution 
must be taken into account in the 0 < 9 < 7t/2 region, to obtain the baryon 
number, B, in the TPM. In fact, when adding the contributions coming 
from the truncated Skyrme model, the Dirac sea and the valence quarks, one 
obtains B = 1 for any value of 9. This can be understood by studying the 
energy of the fundamental eigenstate of the Dirac Hamiltonian [[17]]. When 


9 < 7t/2, the fundamental energy becomes positive and valence quarks leave 
the Dirac sea. In this case, they must be explicitely included. As regards the 
bag energy, the valence quark contribution can be written as 


e q = 3e{9)H(n/2-9), (10) 

where e(9) is the fundamental eigenvalue and H(x) is the Heaviside step 
function. 


Note that our regularization prescription lead to a finite bag energy (de¬ 
pending on Kq and Km it )• All derived quantities will also be finite. Such 
is the case of the axial flux through the boundary of the sphere, 

(j) f (R,9) = [ dQ n^n a (j£ a ) = -^e t (0), (11) 

Jr=R R d9 

where e t {9) = [A e c (9) + e q {6) + Kmit\- The last term, coming from the MIT 
bag, is ^-independent and does not contribute to the axial flux. 
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As it was proposed in reference [T8|], the vanishing of the flux when R —> 
0 determines the Kq renormalized constant. This physical imposition is 
consistent with the analysis we will perform later, when treating the TPM. 


The external phase of the TPM is described by the Skyrme model [||, |j]. 
The lagrangian can be written as 


C = 


^Tr (dpUdPrf) + —~^Tr [(d,U) U\ (d u U) U f 


32e 2 


( 12 ) 


where the scalar field U(x) takes its values in the SU{2) group. 

Two parameters have been introduced here: F 7r , associated with the pion’s 
decay constant (experimental value F% xp = 186 MeV) and e 2 , which repre¬ 
sents the strength of the stabilizing term. These parameters will be adjusted 
later in the framework of the two phase model |2|]. 

The Skyrmion is a topologically stable classical solution of the Lagrangian 
in (0). when the whole space is considered. It is given by 


U 0 (x) = e ie(r)( ™ } , (13) 

where 9(r) is what we are calling the chiral angle. Spatial and isospinorial 
indexes are linked in the argument of the exponential. 

By the imposition of the boundary conditions 


9{r = 0) = 7T 6{r) -> r ^ 0 , (14) 

an skyrmion of topological baryonic number (winding number) B — 1 is 
obtained (in the pure Skyrme model). 

Replacing the Skyrme ansatz in the Lagrangian (|1^), the equation of 
motion is obtained as a nonlinear differential equation for 6{r) 

The Skyrme lagrangian is invariant under SU(2)l ® SU(2)r chiral trans¬ 
formations 

U -> AUB ~ l , (15) 

where A and B belong to SU( 2). When A = B~ 1 , we are in the case of axial 
symmetry leading to the locally conserved axial current. Chiral boundary 
conditions guarantee its conservation even at the boundary, when the TPM 
is considered. 
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The flux of the axial current through an sphere of radius r, in terms of 
the adimensionalized radius r = is given by 



2nF 7r d,9 „ 
e ~dd r 


2 



(16) 


ft is not difficult to show that, when R —> 0 (6 —> it, 9' < oo ) <f) Sk —> 0. 
It is reasonable to extend such behavior to the flux from the inner fermionic 
phase, using this criterium to determine the value of the renormalized con¬ 
stant Kq. 


In the TPM, the Skyrmion is truncated to the exterior of an sphere of 
radius R. In the R —> 0 limit, the pure skyrmion should describe the baryon 
properties. The contribution Kmit/R to the Casimir energy is forbidden 
in such scheme. Then, the validity of the hybrid chiral model, even in the 
R —> 0 limit, imposes Kmit = 0. 


Once our renormalization scheme has been established, a fit can be pro¬ 
posed for the numerically evaluated Casimir energy. Following [18] we pro¬ 
pose 


e c (9)-3X g sin 2 0 = 3|A(| (jf |<«< l 

+C 2 sin 2 9 + C 4 sin 1 9 + C e sin 6 9 + C$ sin 6 9 j , 

where the coefficients take the values 

C 2 = -0.13381 C 4 = 0.05085 

C 6 = -0.01257 C 8 = 0.01241 ' 


(17) 


(18) 


The required vanishing of the axial flux of fermions in the R —> 0 limit 
inmediately leads to 

k q = -c 2 , 

thus eliminating the contribution proportional to (n — 9) 2 for R 0. 

Having fixed the renormalized constants, the energy of the inner phase is 
as shown in Figure |l|. Our results are totally consistent with those presented 
The symmetry of the Casimir energy (dashed line) about 9 = 7t/2 is 


m 
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Figure 1: (Dimensionless) Inner energy of the chiral bag ;-: Casimir 

energy ; — : Total inner energy 

evident. Such property was a priori expectable from simple properties of the 
eigenvalues of the Dirac hamiltonian of the model. 

In the same figure, the solid line represents the total energy of the inner 
phase (when the valence quark contribution is also taken into account). As 
in reference 0, the resulting interior energy is a smooth function. 

In the TPM picture, the axial flux through the boundary of the defect 
should be continuous for all values of the bag radius. So, a fine tunning of 
the strength of the stabilization term e(R) must be performed to ensure 

(,b f (R) =(p Sk (R ) for all R. (19) 

To impose this condition, the knowledge of the chiral angle as a function 
of R , 0(R), is necessary. As it was said, its value can be obtained numerically, 
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by solving a nonlinear differential equation. But, following the proposition 


of M. Atiyah and N. Manton [19], U(x) configurations can be constructed 
by evaluating the holonomy of Yang-Mills fields with topological charge k, 
in the time direction. This is derived from a t’Hooft instanton of width A. 
For the k — 1 case, the resulting chiral angle is 


9(f) = 7T 


1 - 


( 20 ) 


' 1 + (b J 

Replacing 6(f) in the energy of the skyrmion[[7] 

M S k = — J d 3 x Csk 

= r {f {r 2 & 2 + 2s in^) + TfAl { 2rV + *,»«)} . (21) 

and minimizing with respect to A, the value of this parameter is fixed to A = 
1.452 [Hj. The resulting prohle is very close to the one obtained numerically. 
For computational convenience, we will use this aproximation for the rest of 
the paper. 

With the Atiyah-Manton prohle, equation (|l^) gives 


37T 


de t (9) 


e 2 (9) 


32A 2 (vr - 9) 3 


l + 2 6(2n-e) .2 Q 

1 + A 2 (n~e ) 2 &ln U 


do 


( 22 ) 


when parametrized by the chiral angle 6. 

When the 9 —> tt (R —> 0) limit is taken, the stabilization strength for 
the pure skyrmion is obtained 


e(R = 0) = 4.216 , (23) 

to be compared with the value e = 5.45 obtained by G. S. Adkins, C. R. Nappi 
and E. Witten by fitting the masses of the nucleon and the A particle. 

Now, the value of F n can be fixed in the scheme proposed in [0. Hav¬ 
ing e and the nucleon mass M* xp = 938MeV, expression (9) of 0 (with A 
and M calculated for the Atiyah-Manton prohle) leads to F w = 99.59 MeV. 
This value, far from the experimental one, is, however, near F n = 129 MeV 
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Figure 2: Strength of the stabilization term e(R) 

obtained in [[7|. As a consistency check, the A particle mass has been cal¬ 
culated with our parameters, giving Ma = 1206 MeV (experimental value: 
M“ p ' = 1230 MeV). 

The calculations just detailed complete the determination of the param¬ 
eters of the Skyrrne model, in a consistent way with the TPM under study. 
Now, recovering the dimensional variable r, the strength e(R) is as shown in 
Figure |2|. 

To complete our task, the energy contained in the Skyrrne external sector 
must be evaluated. To do this calculation, the integral in (|2T|), with the 
lower limit truncated to the radius of the defect, must be studied. We have 
performed the numerical evaluation of this quantity using the Atiyah-Manton 
profile and the parameters of the Skyrrne model F w and e(R), thus obtaining 
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Figure 3: Two Phase model energy — : Complete model ; - • - : Skyrmion’s 
sector ; — : Bag’s sector 

the energy in the Skyrme sector as a function of the defect radius. 

Figure |3] is the main result of this paper. It shows in dashed lines, from 
top to bottom, the energy of the external Skyrme phase and the energy of 
the defect, as functions of the position of the limit between the phases. Also 
shown, in solid line, is the total energy of the hybrid chiral bag model. This 
last shows a remarkable independence with the bag radius (for 0 < R < 1 fm), 
as suggested by the CCP. 

Summarizing, we have employed the results presented in [11]. [TI| for the 
internal Casimir energies in a chiral bag model. In these papers, by the use 
of analytical regularizations, a renormalized Casimir energy, dependent on 
the Kq and Kmit constants, was obtained. In the present paper, a TPM 
was completed, by introducing an external Skyrme phase. The renormalized 
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constants Kq and Kmit , as well as the parameters F n and e(R ) of the trun¬ 
cated Skyrme Lagrangian were determined according to physical conditions, 
suitable for the TPM. 

In reasonable agreement with the Cheshire Cat hypothesis, the total en¬ 
ergy of this model shows an approximate independence with the bag radius 
(separation limit between the phases) in the range of 0 < R < 1 fm. 


The study of CBM at finite temperature 
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is an interesting effec¬ 
tive approach to the analysis of deconfinement transitions. In those refer¬ 
ences, succesive approximations to the problem, based on the validity of the 
Cheshire Cat hypothesis at T = 0 have been made. The present results give a 
ground to such hypothesis, thus making it sensible to look for the presence of 
deconfinement transitions only in the temperature-dependent contributions 
to the free energy of the bag. 
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